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^- ABSTRACT 

Svi ■ The soft X-ray excess emission in the southern cluster Sersic 159-03 represents 

hiterto the strongest effect of its kind. Emission in the ~ 0.2—0.4 keV passband is 

^ ■ detected far in excess of the expected contribution from the hot phase of the intra- 

^ ■ cluster medium, and extends to the X-ray signal limit of the cluster. Our analysis 

cn ■ of ROSAT PSPC observations reveal that the soft excess can be interpreted as 

Q . either a thermal or non-thermal effect, and the high data quality allows to place 

^—1 . tight constraints on the two currently competing models. However, each model 

->^ . now implies major revisions to our understanding of clusters of galaxies: either 

p . 'warm' gas masses in similar amounts to the hot gas, or relativistic particles in 

Q . or above equipartition with the hot phase, appear to be unavoidable. 

^: 

cd ■ Subject headings: galaxies: clusters: individual (Sersic 159-03) - intergalactic medium 

>: 
^: 

C^ ■ Introduction 

During the five intervening years since the discovery of EUV and soft X-ray excess 
emission in the Virgo and Coma clusters with the Extreme Ultraviolet Explorer (EUVE) 
and i?(95'j4r missions (Lieu et al. 1996a,b), the original results have been cross-checked by 
follow-up observations of the same objects plus a larger sample of nearby galaxy clusters 
(Fabian 1996; Bonamente, Lieu and Mittaz. 2001a,b; Bonamente et al. 2001c; Mittaz, Lieu 
and Lockman 1998; Lieu et al. 1999; Lieu, Bonamente and Mittaz 1999,2000; Berghofer, 
Bowyer and Korpela 2000a,b; Bowyer and Berghofer 1998; Bowyer, Korpela and Berghofer 
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Bowyer, Berghofer and Korpela 1999; Arabadjis and Bregman 1999; Reynolds et al. 
Kaastra et al. 1999; Dixon, Sallmen and Hurwitz 2001; Dixon, Hurwitz and Ferguson 
Valinia et al. 2000; Fusco-Femiano et al. 2000 Buote 2000a,b,2001; Arnaud et al. 



2001), using complementary data from the BeppoSAX, HUT, FUSE and RXTE missions. 



-2- 



Recently, we have undertaken a more extensive study of soft X-ray emission in galaxy clusters 
of low/intermediate N^^ with the aim of assessing the impact and cosmological significance 
of the soft excess phenomenon. Our efforts led to the discovery of soft excess emission in 
a sample of clusters in the Shapley concentration (Bonamente et al. 2001b), and while the 
analysis of more sources is still underway, we report in this Letter the brightest soft X-ray 
(~ 0.2 — 0.4 keV) excess cluster within the ROSAT PSPC database. The excess emission 
reaches ~ 100 % of the hot intra-cluster medium (ICM) contribution, 3-4 times higher than 
the same for any other PSPC targets (e.g., A1795 and Virgo; Bonamente, Lieu and Mittaz 
2001a). The inverse- Compton (IC) interpretation of the soft emission (see section 3) results in 
a non-thermal pressure that greatly exceeds that of the thermal gas: the excess of Sersic 159- 
03 has more demanding energetic requirements than all previously known cases (e.g.. Coma, 
Virgo and A1795; Lieu et al. 1999, Bonamente, Lieu and Mittaz 2001a). Alternatively, if the 
origin is thermal, the mass implications of the putative 'warm' gas are also extreme, with 
the warm-to-hot gas mass ratio ~ 50 % (this figure is matched only by the emission seen 
by EUVE for A2199, however at lower statistical significance; Lieu, Bonamente and Mittaz 
2000). 



X-ray, infrared and 21 centimeter observations of Sersic 159-03 

Sersic 159-03 is a rich southern galaxy cluster, also known as Abell SI 101, at a redshift of 
z=0.056 (Abell, Corwin and Olowin 1989). The cluster's X-ray luminosity {Lx = 5.35 x 10'^^ 
erg/s in 0.5-2.0 keV, De Grandi et al. 1999) is typical for a cluster of its temperature (e.g., 
Wu, Xue and Fang 1999), and the emission has a radial extent of ~ 10 arcmin (Fig. 1). 
Recent XMM observations (Kaastra et al. 2001) show a 'cooling' of the hot ICM in the 
central regions and no obvious peculiarities in the X-ray morphology, consistent with our 
PSPC data. The cluster, moreover, has no known radio emissions associated with it. 

By means of three pointed PSPC observations of this source, taken between May 1992 
and May 1993 with a total exposure of about 19,000 seconds, we found a very bright excess 
emission in PSPC's 1/4 keV band (also known as R2 band, ~ 0.2 — 0.4 keV). Fig. 2 shows 
how detected emission (light blue, dark blue and black crosses) compare with the expected 
emission from the hot intra-cluster medium (ICM, solid lines), e.g. as described by the recent 
X-ray analysis of Kaastra et al. (2001). The three observations were modeled simultaneously, 
in order to reduce statistical uncertainties. 

As extragalactic soft X-ray emissions are absorbed by the interstellar medium, an accu- 
rate measurement of the line-of-sight Galactic hydrogen (HI) column density {Nh) is crucial. 
The Dickey and Lockman (1990) HI maps at resolutions of respectively 1, 4 and 9 square de- 
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grees centered at the cluster's position indicate a region of sky free from spatial HI gradients, 
with an HI density at the cluster position of 1.79 xlO^*^ cm~^. Since Galactic infrared 100 /im 
emission correlates well with HI 21-cm radiation (Boulanger and Perault 1988), we further 
consulted IRAS 100 /im maps (Wheelock et al. 1994) at a finer resolution of 3 arcmin. The 
100 /im emission is very smooth over the cluster's region, and it confirmed our Nh value 
for the cluster's Galactic HI. This value is henceforth used in the spectral analysis ^. We 
refer to a previous paper (Bonamente et al. 2001b) for the details on PSPC data analysis 
techniques. 

Spectra of azimuthally averaged annular regions were modelled with a single tempera- 
ture optically-thin plasma emission. ^ The Galactic absorption was modelled with the codes 
of Morrison and MacCammon 1983 (MM83; 'WARS' in XSPEC language) and of Wilms, 
Allen and McCray 2000 (WAMOO; 'TBABS' in XSPEC). The He cross-sections of MM83 are 
in good agreement with the recent compilation by Yan, Sadeghpour and Dalgarno (1998), 
as also shown in a previous paper (Bonamente et al. 2001a). The WAMOO code provides 
however the most up-to-date cross-section compilation for absorption of X-rays in the in- 
terstellar medium. In the energy range of interest, 0.2-0.4 keV, the two codes are nearly 
indistinguishable at the resolution of the PSPC (see Fig. 4 of WAMOO), and yielded sta- 
tistically consistent results, see Tables 1 and 2. The emission code employed was MEKAL 
(after Mewe, Gronenschild and van den Oord, 1985; Mewe, Lemen and van den Oord 1986; 
Kaastra 1992). 

When applied to the entire passband of the PSPC (0.2-2.0 keV), the fit is formally 
unacceptable (e.g., red crosses in Fig. 3), and a region of depleted flux (0.5-1.0 keV) appears 
to accompany the 0.2-0.4 keV excess. The fit is satisfactory when applyed only to energies > 
0.5 keV, with the soft excess then abundantly revealed at E< 0.4 keV (green crosses. Fig. 3). 
The best-fit parameters of the hot phase agree with the most recent X-ray XMM observation 
(Kaastra et al. 2001). In Table 1 we summarise the results. We emphasize that the soft 
excess emission is so strong that, in order to be interpreted as the low-energy tail of the 
hot ICM emission, a Galactic column density as low as 5-8 xlO^^ cm~^ would be required, 
which is not only severely discrepant from the measured values, but even lower that that of 
the 'Lockman hole' where N^ reaches global minimum. The soft emission cannot therefore 
be attributed to peculiarities in the Galactic HI distribution. 



^The low declination (^ —42°) renders this source not easily observable with the NRAO 43-mt Green 
Bank telescope we employed in our previous analyses. 

^The spectral analysis was performed with the XSPEC software. PSPC pulse-invariant (PI) channels 
1-19 and 202-256 are not used, as potentially affected by calibration uncertainties. 
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In Figure 4 we show the trend of the PSPC excess. The emission extends to a radius 
of 800 kiloparsec ^ (=9 arcmin); in the 9-12 arcmin region, where only marginal cluster 
emission is detected (Fig. 1), the data point is consistent with zero flux, thereby confirming 
the integrity of our background subtraction technique (Bonamente et al. 2001b for reference). 
As comparison, we note that from the XMM/EPIC data of Kaastra et al. (2001) evidence 
for X-ray emission is apparent out to a radial distance of about 10 arcmin, in agreement with 
our PSPC data. The soft excess emission therefore persists to the X-ray detection limit. To 
illustrate the significance of the role played by the model uncertanties in our determination 
of the 1/4 keV contribution from the hot ICM, we consider as example the 3-6' region. When 
the statistical error in the 0.5-2.0 keV PSPC data (the portion of the spectra used for model 
determination) is taken into account, the error in the fractional eccess increases only from 9 % 
(see Fig. 4) to 10%. Similar conclusions apply to the other regions, indicating that evidence 
for a soft X-ray component from Sersic 159-03 is unaffected by statistical uncertanties in the 
hot ICM modelling ^. 



Interpretation 

How do we interpret such a strong signal? A 'warm', optically-thin phase of the intra- 
cluster medium (ICM) at sub-megakelvin temperatures can in principle accont for the excess 
emission. Thermal conduction can be suppressed, either for the presence of tangled magnetic 
fields (Chandran et al. 1999) or if the warm phase is sufficiently 'clumpy' (Cowie and McKee 
1977; Smith and Lilliequist 1979), yet the time scale for radiative losses of a 'warm' gas is 
sufficiently short to necessitate rapid and continuous replenishment (such as that of the 
'mixing layers' scenario of Fabian 1997), if the warm phase existed in a steady state for 
cosmological times. In the neighborhood of T ~ 10® K, the cooling time of a gas with solar 
abundances is tcooi ~ 3 x 10^ x nZl years, where n_3 is gas density in units of 10~^ cm~^ 
(Landini and Monsignori-Fossi 1990), the characteristic density of the warm phase in the 1-3' 
and 3-6' regions for a volume filling factor / = 1(0</<1). Our observations however 
determine the so called emission measure of the warm phase {EM oc J ri^dV, V being the 



■^A Hubble constant of Hq = 50 km s ^ Mpc ^ is assumed here and after. 

^Krick, Arabadjis and Brcgman (2000) argue that sub-arcmin variations of N^ could in principle account 
for a 20-30 % excess. Although the relevance of this phenomenon to our Sersic 159-03 observations is 
uncertain, the strength of the PSPC excess here reported largely exceeds that figure. On occasions where 
comparison with stellar Lyman-a and quasar X-ray spectra could be made, the agreement among the different 
methods of determining Nh would suggest an error less than ^ 10^^ cm~^ (see Lieu et al. 1996b and 
references therein). 
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voluine of the emitting gas), and / remains unknown. For a given detected EM, the total 
warm gas mass M in a given region is a function of /, M oc ^/f; as the density increases 
with decreasing / (n oc 1/y/f), the coohng time accordingly decreases as tcooi oc v7- The 
mass implications are very demanding: this new component will have a mass comparable 
to that of the hot ICM, see Table 2, and those budgets can be alleviated if / is very small, 
however at the price of a correspondingly shorter cooling time. 

Alternatively, an inverse Compton (IC) origin of the emission can be advocated (Sarazin 
and Lieu 1998). Diffusive shock acceleration (Axford et al. 1977; Bell 1978a,b; Blandford and 
Ostriker 1978) may be responsible for a population of relativistic electrons with a power-law 
spectrum in the ICM, and its IC mechanism preserves a similar power-law shape, whereby 
the two differential number distributions of electrons (N(E)oc E^'^) and emitted photons 
(L(e) oc e~") are related by the equation a = (1 +/i)/2. Relativistic electrons of 7 =300-700, 
required to emit EUV or soft X-ray through IC scattering, straddle spectrally between low- 
energy Coulomb losses and high-energy radiative losses (Lieu et al. 1999; Sarazin 1999), and 
may survive for a significant fraction of a cluster's lifetime. Here we fit the photon index a to 
the data. Table 2, with the exception of the 3-6 arcmin region where a it is fixed at a=1.75, 
the expected value from a cosmic ray (CR) population of fi=2.5 in accordance to Galactic 
CR. The modelling is formally acceptable, although with a spectral index somewhat steeper 
than that of Galactic CR. However, serious problems confront the energetic requirements of 
the IC model: the pressure of CR electrons with 7 in the range 475-700 (those which emit IC 
radiation in the E=0.2-0.4 keV band) exceeds that of the hot ICM by a factor of order 2-3! 
Moreover, this pressure estimate is necessarily a lower limit, as higher- and lower-energy 
electrons may be present, as well as CR ions which, although not visible, will further drive 
the CR budget to absurde values (Lieu et al. 1999; see also Miniati et al. 2001). 

As often happens in rich clusters of galaxies, the X-ray centroid, here coincident with 
the central cD galaxy ESO 291-9, is associated with a lower temperature for the hot ICM 
(Kaastra et al. 2001). Such a temperature decrease is also clearly visible in our PSPC data. 
Yet the detected soft X-ray excess cannot be directly related to it: firstly, because the excess 
covers an area about 25 times larger than the region affected by this 'cooling' (which has 
a radius of about 1.8 arcmin, Kaastra et al. 2001; Allen and Fabian 1997); secondly, and 
more importantly, because we already accounted for the effect in our data modeling (see 
best -fit temperatures in Table 1). The behaviour here is therefore analogous to, e.g., that 
of the Virgo and A1795 clusters, where the excess emission spreads over a much larger area 
than that of the central cooler region (Lieu et al. 1996a; Bonamente et al. 2001a). The 
central galaxy of Sersic 159-03, which also contains an IR source, is well contained within a 
radius of 1 arcmin (Hansen et al. 2000): the soft excess emission must then be a genuinely 
cluster-wide phenomenon. 
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Discussion and conclusions 

The high statistical significance of the soft excess emission of this cluster confounds both 
thermal and non-thermal models, Table 2. Either model can fit the data, and calls for an 
overhaul in our understanding of galaxy clusters. 

Following the non-thermal interpretation, the large pressure budgets of Table 2 cannot 
easily be mitigated; on the contrary, they are strict lower limits to the relativistic particle con- 
tent in the cluster. The presence of a population of relativistic particles above equipartition 
with the hot gas clearly has implications on a cluster's hydrostatic balance (e.g., Berezinsky 
et al. 1997), and the role of intra-cluster CR on the cluster's evolution (e.g., heating of the 
ICM gas) could be much more significant than previously thought. At present it is not known 
whether strong shocks, e.g., those induced by cluster mergers, may accelerate such a large 
amount of CRs. It is possible that the cluster experiences many such shocks throughout its 
lifetime; the ICM would then accumulate those CR electrons with the longest radiative life 
time, precisely the 7 ~ 300-500 electrons which emit in EUV and soft X-rays through IC 
scattering (Sarazin and Lieu 1998; Lieu et al. 1999). 

On the other hand, large warm gas masses could be sustained at the interface of cold 
gas clouds and the hot ICM (Fabian 1997), the presence of a 'cold' (T < 10^ K) phase in 
the ICM is in fact an interpretation of the PSPC images of the Coma and Virgo clusters 
(Bonamente et al. 2001c). In the absence of such a replenishment mechanism, however, the 
very short radiative cooling time renders the thermal model untenable. 

The bright soft excess of Sersic 159-03 here reported presents a new and unavoidable 
reality: whether the ultimate explanation be thermal, non-thermal or some other origin, a 
major effect at work in the intergalactic medium has hiterto been completely ignored. 
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Table 1: Spectral modelling of PSPC spectra of Sersic 159-03 with the simple 1-T model. 
Here and after errors are 90 % confidence (x^ + 2.7 method), where no errors are reported 
parameters were held fixed, according to the X-ray study of Kaastra et al. (2001). Fit 
is performed by minimizing the x^ function, d.o.f. is degrees of freedom. 'A' is elemental 
abundance following the relative proportions of Anders and Grevesse (1989). Photoelectric 
absorption code is MM83, in square brackets the new WAMOO code is employed, with same 
results as MM83. 



Region 




0.2-2.0 kcV fit 






0.2-2.0 keV 


fit. frco NH 






0.5-2.0 kcV fit 




(arcm.) 


kT 
(kcV) 


A 


X^ (d.o.f) 
(koV) 


kT 


A 
10™ cm-2 


NH 


X^(d.o.f) 


kT 
(koV) 


A 


X^ (d.o.f) 


0-1 


1-44 ±o:o2 


0.08 ±8:81 


264(167) 


1.88 ±p? 
1.92 ±l\l 


n / f 0.12 

0-4 ±0.09 


0.95 ±°il 

0.92 ±8i| 

0.52 ±o;i| 

0.51 ±8:« 


143(166) 


1.87 ±«:lf 
1.87 ±8:11 


0.41 ±«:;^ 

0.41 ±8:}* 


97(137) 




[ 1.48 ±0.07 


0.08 ±8:81^ 


272(167) 


0.43 ±8}'' 


142(166) 


97(137) 1 


1-3 


1.38 ±0.08 

[ 1.38 ±«:f 

[2°t^4 

I 2 ±0.76 


0.03 ±8;8i 
±8°= 


335(162) 
342(162) 


3 ±0.5 


n B 1 0.34 

05t*-%^ 
0-0 ±0.21 

0.6±8.5 


154(161) 
153(161) 


2-9±8:? 

2.9 ±0.65 


5±''-37 

0-J ±0.23 
0.6±8.5 


115(132) 
115(132) ] 


3-6 


0.02 ±0:81 


111(128) 


4.o±?:| 


o-7±P„ 


88(127) 


3.8±f.4 


66(99) 




±8-"i 


112(128) 


4.1 ±?;1 


0.77 ±8.6 


0.68 ±8ii 


88(127) 


3.7 ±M 


0.64 ±8.63 


66(99) ] 


6-9 


2 


0.3 


67(84) 


2 


0.3 


l-3±0.65 


65(83) 


2 


0.3 


31(55) 




[2 


0.3 


67(84) 


2 


0.3 


1 ., 1 0.85 
1.3 ±0.4 


66(83) 


2 


0.3 


31(55) 1 
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Table 2: Spectral modelling with a 2-T model ('hot' and 'warm' phases of ICM) and with a 
non-thermal model. Mass implications for the former and pressure budgets for the latter are 
reported; / is filling factor of the warm component. See Table 1 for explanation of Galactic 
absorption codes employed. 



Region 
(arcm.) 



(kcV) 



2-T model 
kT„ 
(koV) 



X^(d.o.f) 



M^/M^ 



0-1 


[ 1-9 ±S:f2 
2-9 ±8:Ii 


0.43±i5;E? 


1-3 


0.64±00^P 
0.53±g:tl 




I3±8:f 


0.67 ±0-^2 


3-6 


3.5±ll 


0.3 




[3.5±?i 


0.3 



0.078±lJ;8?i 
0.07±0.035 
0.06±8:»ti 
0.07 ±[]:0f 
0.06±0.0035 



130(165) 
129(165) 
142(160) 
141(160) 
88(127) 
88(127) 



V7x 0.38±J];gg 
y/Jx 0.34 ±0-fe 
VTx 0.55±p 
^x 0.33 ±Jj-|| 
v7x 0.46±J5 
v7x 0.61 ±|!i^42 



kT 
(koV) 



Non-thermal model 

X^(d.o.f) P, 



t/ier 

10" 



Pic 

(erg cm ) 



■94±8:ii 
28^8:^' 

2.8±n.4 
2.8±0:| 


0.5 


., 91 0.9 
3.2±j,2 


136(166) 


0.5 


3.2±H 
4.8±S;i 


135(166) 


0.4 


148(161) 


0.4 


4.7±J:i 


148(161) 


2.5 


0.2 


1.75 


96(129) 


2.5 


0.2 


1.75 


96(129) 



7 
0.7 
0.14 



22±i3 
20 ±11 ] 

2.3 ±I:i 

2 ±1,3 1 
0.33 ±0;h 
0.36 ±0.12 1 
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Figure 1: Radial profiles of X-ray brightness in PSPC R2 (0.2-0.4 keV, red) and R37 
(0.4-2.0 keV, green) bands in units of counts/s/sq.arcmin, the lines mark the corresponding 
background levels. Background is determined from an annulus ~ 40' off-axis. 

Figure 2: Three PSPC spectra of the central region (radius 3 arcmin) of Sersic 159-03. 
The data, available through the public archive (identification numbers 800397, 600182 and 
600429al), are individually rebinned to obtain a minimum S/N > 4 in each spectral channel. 
The solid lines describe the best-fit single temperature model (see text for description) 
that highlights the low-energy excess of photons. Agreement among the 3 observations is 
excellent, see bottom panel. 

Figure 3: Coadded spectrum of 1-3 arcmin region of the cluster (black crosses) overlaid 
on the best fit 1-T model obtained by fitting the whole PSPC band (left, red solid line) and 
by fitting only energies > 0.5 keV (right, green solid line). Residuals (red and green crosses) 
reveal in both cases excess of soft photons at energies < 0.4 keV. 

Figure 4: Fractional excess in 1/4 keV PSPC band (0.2-0.4 keV) as function of radial 
distance from the cluster's center. Fractional excess is defined a.s rj = {m — p)/p, where m is 
measured soft flux (1/4 keV band) and p its prediction according to the 0.5-2.0 keV single 
temperature model. The 9-12' region was modelled with the same parameters of the 6-9' 
annulus (Table 1); vertical semi-diameters are l-cr statistical errors. The excess component 
ranges between 50 and 100 % of the hot ICM gas in the soft band, the strongest observed 
to date in PSPC data. 
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